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Abstract 

The dynamics of a soft sphere model glass, studied by molecular dynamics, is investigated. 
The vibrational density of states divided by uj 2 shows a pronounced boson peak. Its shape is in 
agreement with the universal form derived for soft oscillators interacting with sound waves. The 
excess vibrations forming the boson peak have mainly transverse character. From the dynamic 
structure factor in the Brillouin regime pseudo dispersion curves are calculated. Whereas the 
longitudinal phonons are well defined up to the pseudo zone boundary the transverse ones rapidly 
get over-damped and go through the Ioffe-Regel limit near the boson peak frequency. In the high 
q regime constant-a; scans of the dynamic structure factor for frequencies around the boson peak 
are clearly distinct from those for zone boundary frequencies. Above the Brillouin regime, the 
scans for the low frequency modes follow closely the static structure factor. This still holds after a 
deconvolution of the exact harmonic eigenmodes into local and extended modes. Also the structure 
factor for local relaxations at finite temperatures resembles the static one. This semblance between 
the structure factors mirrors the collective motion of chain like structures in both low frequency 
vibrations and atomic hopping processes, observed in the earlier investigations. 

PACS numbers: PACS number(s): 61.43.Fs, 63.50,+x, 66.30.Fq 
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I. INTRODUCTION 



Despite large efforts, the vibrational dynamics of amorphous materials and glasses is 
still not fully understood and subject to a controversial debate. At a first sight there are 
many similarities between glasses and amorphous materials and their ordered counterparts, 
the crystals. The vibrational densities of states are normally similar and typical average 
frequencies or Debye temperatures do not change dramatically from one state to another. 
This reflects the the similarity of the short range order of both states, primarily the average 
nearest neighbour distance is similar. At higher, but still low, frequencies the vibrational 
density of states (DOS), Z(u), in glasses exceeds the Debye value and a maximum is seen 
in the inelastic scattering intensity, the so called boson peak (BP). Due to disorder the 
sound waves are increasingly damped with increasing frequency and eventually the damp- 
ing exceeds the Ioffe-Regel limit: The mean free path of the phonons diminishes below 
their wavelengths, the phonons are over-damped. The vibrations can still be understood 
by harmonic theory but their eigenvectors will have a complicated structure reflecting the 
structural disorder. The frequency, where the Ioffe-Regel limit is reached, is normally 
near the BP frequency, uj-qv- As in defect crystals, one finds in addition to the "band" modes 
localized vibrational modes above the maximum frequency of the band modes or in gaps of 
the spectrum. Additionally, particularly in molecular glasses, there can be the equivalent of 
the optical vibrational modes of crystals, albeit often broken up by disorder. 

The dispute centers mainly on two issues: first, the nature of the excess of low frequency 
modes above the Debye value, i.e. the origin of the BP, and secondly on the properties of the 
phonons at higher frequencies. The term phonon is used here for an extended vibrational 
mode which can be classified by a g-value. The latter question has gained great interest 
by the advance of modern synchrotron sources which allow to measure structure factors in 
Brillouin scattering. 

The frequency domain, typically around v — 1 THz, of the BP is accessible to a variety 
of spectroscopic methods like Raman [1] and neutron scattering |2j. From the temperature 
dependence of the measured intensities it has been deduced that the excitations are essen- 
tially harmonic vibrations. The boson peak is, therefore, a maximum in Z(uj)/uj 2 . Mostly 
it is not a maximum of the vibrational spectral density, Z(u), itself. 

Over the years a large number of explanations of the excess intensity and the BP have 



2 



been proposed. First, there are some materials where some molecular unit has an inherent 
low frequency mode, such as in plastic crystals. These modes will trivially show as "excess" 
intensity in Z(u). We will not consider these. Secondly, there are materials where the BP is 
a real effect of disorder effect. The model, studied here, belongs to this class it has further 
no molecular vibrations or low lying optic modes in its crystalline form. Some authors 
attribute the BP to vibrations of clusters of atoms of typical sizes . The origin of these 
clusters is unclear and they have not been identified in numerical simulations. Another 
popular qualitative explanation of the Boson peak is a softening of acoustic phonons by 
static disorder 0, |f| due to elastic Rayleigh scattering. However, estimates show that the 
Rayleigh contribution is too small to explain the experimental data on thermal conductivity 
in glasses (|. Sometimes the BP is related to low lying optic modes of parental crystals p, 
y,0]. Whereas the BP is a general feature of glasses such crystal structures with soft optic 
modes cannot be identified always. Such a mechanism is possible in some cases. A distinction 
optic modes, which disorder has broken into segments, and resonant vibrations is not very 
clear in all cases. 

Recent work on harmonic lattice models demonstrated that softening of disordered force 



constants can smear and push to low frequencies peaks which exist in the crystalline DOS 
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Q|. In another approach the vibrations of a random distribution of atoms, interacting 
with a Gaussian-shaped pair potential, was studied in a harmonic scalar approximation. 

A different explanation is offered by the soft potential model (SPM) Q|- This model 
gives a unified description of the glassy low energy dynamics: tunneling, local relaxation 
and excess in low frequency vibrations. For the latter it postulates low frequency oscillators 
which interact bilinearly with the sound waves and hybridize with them. Thus quasi-local 
vibrations, also called resonant vibrations, are created jl^l ]. Fitting the model to the ex- 
perimental data, one finds effective masses of 20-100 atomic masses for the entities moving 
in these effective soft potentials Q.The sound waves in turn will be damped by resonant 
scattering with these modes [3]. For low frequencies the SPM predicts an increase of the 
excess of the DOS as Z excess (uj) oc u A . This rapid increase of Z excesB (uj) will lead to level 
repulsion and an eventual cross over to Z excess (cc;) oc uj [l9|. Recently it has been shown 
that the interaction between the soft oscillators in harmonic approximation causes unstable 
vibrations, negative eigenvalues. The structural configuration becomes unstable. The re- 
sponse of the glass is then a local relaxation to a new stable configuration which always exist 
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due to anharmonicity. The harmonic spectrum in this new stable configuration exhibits a 
BP with a universal shape as it is observed in experiment. In particular the shape of the 
BP is independent of the actual magnitude of the anharmonicity j^J. 
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22 1 confirmed the existence of low 



Computer simulations of a soft sphere glass (SSG) 
frequency quasi-localized vibrations with effective masses ranging from ten atomic masses 
upwards. The modes were centered at structural irregularities with large local strains j^j]. 
Regions of local strain have also been observed in earlier computer simulations l24fl . Similar 
effective masses have since been observed in many different simulations, e.g. Si0 2 [25j], Se 26 1 



in Ni-Zr 1271, Pd-Si \2i 



Ni-B [^S I, amorphous ice and amorphous and quasi- crystalline 
Al-Zn-Mg |3l(. In an earlier simulation of amorphous silicon low frequency localized vibra- 
tions have been observed at coordination defects 3^|. It should be noted that computer 
simulations with the usual periodic boundary conditions have a lower cutoff, ^min-sound for 
the sound waves, which is given by the maximal wavelength fitting into the periodicity 
volume. Quasi localized vibrations with frequencies uj < u; m m_ S ound appear, therefore, in 
simulations as truly localized vibrations which can be used as an easy way to identify them. 
For uj > cumin-sound hybridisation with the sound waves occurs and the local vibrations are 
delocalized, they become quasi-localized. This will lead to a small shift in Z exceSB (u) but 
does not affect the number of excess modes. For a more detailed discussion see [3c 
For frequencies uj > ^Vm-sound it is still possible to change from the basis of exact harmonic 
eigenvectors to a basis of sound waves and local vibrations with a small bilinear interaction 
|34j |. The latter basis is the one used in the SPM. 

High resolution synchrotron scattering made it possible to observe the dynamic structure 
factor, S(q, uj), in Brillouin scattering 35]. For a recent review see j^. As in crystals, At the 
lowest frequencies glasses have well defined longitudinal and transverse sound waves. Only 
the longitudinal one is observed in Brillouin scattering at low q. At the lowest frequencies the 
inelastic scattering function or the dynamic structure factor shows, as function of frequency 
for a given q, only one peak which is ascribed to the longitudinal sound wave. At higher 
frequencies a second peak, ascribed to the transverse sound is observed j^J]. This can be 
understood considering that in a disordered system a transverse sound wave acquires a 
longitudinal component and vice versa. Fitting a damped harmonic oscillator to S(q, uj) in 
amorphous silica 



I and glycerol [39|] , one group concluded in that propagating collective 
excitations exist in glasses at high g-values, i. e. frequencies far above the BP. From 



4 



measurements of densified silica another group 0, III reach the opposite conclusion that 
the BP marks the end of the propagating acoustic modes. 

This question can at present not be answered unambiguously from experiment. Addi- 
tional insight can be gained by computer simulation. Molecular dynamics simulations of 
amorphous silica gave structure factors in good agreement with experiment. A fit with a 
damped harmonic oscillator function gave a damping of the longitudinal mode 7(g) oc q 2 for 
frequencies above the BP (4^ Q . 

The aim of the present paper is to use a simple model glass, representative for dense 
packed metallic glasses, to see what information can be gained from the dynamic structure 
both in the Brillouin regime and for higher g-values. For this we extend our previous study 
of the soft sphere glass 

HQS- 

This glass has the advantage that there is no parent 
crystal with optical modes and no torsional motion as one has in silica which is built from 
more or less rigid Si04-tetrahedra. In the earlier studies j^J we found in a normal mode 
analysis that the Ioffe-Regel limit is reached for the transverse acoustic modes near the BP. 
For the lowest frequency in the calculation the width of the phonon line could be split into 
two approximately equal parts attributed to static disorder and to resonance scattering with 
QLV. 



II. COMPUTATIONAL DETAILS 



The soft sphere glass (SSG) is described by an inverse sixth-power potential 

u{r) = e(^)\A^-) + B. (1) 

To simplify the computer simulation and normal-mode analysis, the potential is cut off at 
r/a = 3.0, and then shifted by a polynomial, A(r/a) A + B, where A = 2.54 ■ lCT 5 e and 
B = —3.43 • 10 _3 e were chosen so that the potential and the force are zero at the cutoff. 
This form of the shifting function was chosen so that its effect is negligible near r/a = 1.0. 
Quantities such as the pressure and average potential energy will be changed by a few percent 
as a result of this truncation, but any changes in the equilibrium structure will be small. 
Without loss of generality one can set e = a — m — 1. Where we do not explicitly state the 
units we infer these "system units" . Note that this soft sphere potential is much softer than 
the often used one with u 12 oc (a/r) 1 2. 
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The calculations were done using the samples prepared earlier for investigations of vibra- 



tions 



and relaxations 



44]. These had been prepared by the following procedure. Liquid 
configurations of 5488 soft spheres were produced via constant-energy molecular dynamics 
(MD) simulation with cubic periodic boundary conditions at a density, pa 3 = 1.0, and tem- 
perature, kT/e ~ 0.54 (about 2.5 times the melting temperature at this density j^j]). For 
the simulation, we used the velocity- Verlet algorithm with a time step of 0.04 - in units of 
(mcr 2 /e) 1 / 2 . The liquid was first quenched within the MD simulation by velocity rescaling to 
a reduced temperature of about 0.005T r The quench rate was about 0.015k/ (mcr^e) 1 ' 2 . Af- 
ter the MD quench, each sample was heated to 0.05T g and aged for several 1000 further MD 
time steps to stabilize the potential energy and to avoid spurious minima. Each system was 
finally quenched to zero temperature using a combination the steepest-descent /conjugate- 
gradient algorithm of our version of the Harwell code DEVIL, the metals version of the 
HADES code. To improve statistics 11 different configurations of 5488 atoms, created in 
this way, were analyzed. 

III. VIBRATIONS 

The vibrational dynamics was calculated in classical harmonic approximation from the 
force constant matrix of the T = minimum configuration. The numerically exact mini- 
mization of the potential energy prevents the occurrence of spurious unstable modes. The 
elements of the force constant matrix are given by 

0M|R"-R"|) 

Different from our previous work all eigenvalues, (cu cr ) 2 , and eigenvectors, e n (a), were 
calculated by direct diagonalisation. 

To localize the BP we calculate the vibrational density of states from the frequencies of 
the 3N — 3 vibrational modes a as 

ZM = (3^X>("-u0) (3) 

where 5 is the discretized S function and () stands for the averaging over configurations. 
Within the numerical accuracy the resulting spectrum is unchanged from Fig. 1 of j^| . The 
maximal frequency is o; max ~ 12.6. There is a slight dip around uo = 9.5 which is charac- 
teristic for the chosen soft sphere glass. It reflects the factor three between the longitudinal 
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FIG. 1: Vibrational spectrum divided by oj 2 versus co. The solid line gives the simulation results, 
corrected by a Debye contribution at the lowest frequencies. The dashed line is a fit with the model 
of Ref . 

and transverse sound velocities. It disappears for harder models, i. e. a larger exponent of 
the cx/r-term in Eq. ^ 

Fig. d (solid line) shows the calculated Z(lo)/lj 2 , corresponding to the BP. From the size 
of our simulation cell we calculate from the elastic constants w m in-sound = 0.52 and 1.76 for 
the longitudinal and transverse sound, respectively. The simulated density of states has no 
sound waves below these cutoffs. To compensate for this we added for the lower frequencies 
a Debye contribution to Z(u), calculated from the elastic constants. The dashed line shows 

n 

the theoretical estimate for the model of interacting soft oscillators and sound waves [20] 



Zthcor(^) — Zj)(u) + Z e 
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where Zn(u) is the Debye spectrum and Z cyLCCSS (uj) the excess causing the BP. There is very 
good agreement between the two curves up to to = 2ubp- The deviation at higher frequencies 
has to be expected. We arrive, in agreement with the earlier work, at an estimate lvbp = 0.5 
which is near to ct'min-sound of the transverse branch and well below the corresponding value of 
the longitudinal branch. The considered soft sphere glass, therefore, represents an interesting 
example of a material where one has a low lying pronounced BP and a strong separation of 
the two branches of sound waves. This distinguishes it from the Lennard- Jones glass. 
In classical harmonic approximation the dynamic structure factor can be written as 



S{q,u) 
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£(q-e»)exp(*q-ir 



6{lo -co a )). 



(5) 



where 6 is the discretized S function and () stands for the averaging over configurations 
and directions of q. Here we omitted the temperature dependent factor kT/hu, the high 
temperature limit of the phonon occupation number. The inelastic scattering intensity is 
given by 

kT 

I(q,u)(x — S(q,u). (6) 

UJ 

Often the structure factor is defined including the temperature factor, it is then directly 
proportional to the inelastic scattering intensity. In our previous paper (^J we used a 
definition excluding the factor 1/u in Eq. 

S FT (q,uj) = uS(q,u). (7) 

This corresponds to a projection of the Fourier transform of the eigenvectors. For broad 
distributions in uj the different definitions will lead to shifts in the maximum and with. 

In analogy with the usual structure factor which projects out the longitudinal components 
of the vibrations one can define a transverse one as 

2 \ 



^(qx e"( ( x))exp(iq-R n ) 



5(u-cj a )). (8) 



It gives, for small q, access to the transverse vibrations which are not observable in the 
scattering experiments. We study the structure factors both in the Brillouin regime, i.e. 
below the first peak of the static structure factor (in our system at gpp = 7cr _1 see Fig. |SJ), 
and above for large q values. 



A. Brillouin regime 

As done by other authors we do constant-g scans of the dynamic structure factor and 
define the maxima in uj as phonon frequencies and derive this way a pseudo dispersion curve. 
This is shown in Fig. |21 for both the longitudinal and transverse branches determined from 
Eqs. |S] and |Hl respectively. The error bars indicate the full width half maximum. 

One sees a very different behaviour of the two branches. Whereas the longitudinal one 
is well defined well beyond the pseudo zone boundary gzB = 0.5 • gpp, the transverse one 
becomes rapidly over-damped. The different broadenings of the two branches indicate a 
large difference in their coupling to the BP-modes. This is seen even more clearly when 
we define the phonons from the scattering intensity, i.e. include the factor kT/u. The 
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FIG. 2: Pseudo dispersion curves in the soft sphere glass (longitudinahtop transverse:bottom) 
calculated from the structure factor without the phonon occupation factor. The error bars indicate 
the full width at half height. The dashed line gives the linear dispersion given by the sound 
velocities, calculated from the elastic constants of the glass. 

longitudinal branch is only slightly affected whereas the transverse one broadens even more 
and above q = 1.8<t _1 the maximal intensity is no longer at some ill denned phonon line but 
drops to the BP frequency. This indicates a transverse nature of the BP. It should be noted 
that the frequencies of both branches at the pseudo zone boundary are much higher than the 
BP-frequency c^bp- This certainly does not justify an explanation of the BP in terms of soft 
zone edge phonons. For low q, both branches show a linear dispersion according to oj — cq. 
The transverse branch shows above q = 1, bo^ 1 an apparent slightly positive dispersion as 
has been seen for the longitudinal branch in SiC>2 42]. 

Fig. El gives the broadening of the pseudo dispersion for different weightings with l/u>, 
corresponding to the definitions via the scattering intensity (Eq. EJ), the correlation function 
(Eq. EJ) or a Fourier transform of the eigenstates (Eq. Ej). As to be expected the definition 
does not change the broadening at the low g-values but is essential above the Ioffe-Regel 




FIG. 3: Halfwidth-halfmaximum divided by lj versus omega, calculated from the dynamic structure 
factor with different w-weighting. Open and full symbols longitudinal and transverse branches, 
respectively; squares: weighting according Eq. [5J circles: Eq. EJ triangles: Eq. [7| The dashed line 
indicates the Ioffe-Regel limit T(q) = uj(q)/2ir. 

limit when T(u) > uj/2tt. The longitudinal branch goes through cufr at q — 4<t _1 beyond 
which value it rises rapidly as q — > qpp. The transverse branch goes through u;ir already 
at q — 0.5a -1 . There are three apparent regions in T(q), best seen in the values derived 
Sft(<7, <^)- At the lowest q the width, T(q), increases with a power of q 2 or larger, followed 
by r(g) oc q oc uj(q) and finally there is an upturn for q — > qpp. 

The soft potential model predicts T(q) oc q 4 for u(q) < u>bp which is derived from the 
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Unfortunately, the low 



initial increase of the excess vibrational DOS Z cxccss (uj) oc u 4 
value of cjbp hi the investigated soft sphere glass excludes this regime from our simulation. 
The simulation cell would have to be large enough not only to allow small g-values it would 
also have to accommodate several mean free paths of these phonons. Otherwise the phonons 
would be scattered by a periodic repetition of equal defects instead of a distribution of 
different defects. For q < la^ 1 we observe for the transverse branch an increase T(g) oc q 2 as 
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has been observed earlier for strongly frustrated systems [46] . The section where T{q) oc ujq 
is broadly compatible with the increase of the total Z{uj) oc uj for uj > ujbp- Above the BP 
the transverse modes are strongly mixed with the quasi-localized vibrations and amongst 
each other. Therefore all transverse modes can be considered as scattering centers. 

Our results for the pseudo dispersion curves agree with the earlier calculation by Caprion 



et al. [47[. There is, however, a marked difference in the calculated widths. This might be 
due both to the better resolution in the present work and also to the different procedure used 
to calculate T{uj). We define the width directly from the width at half maximum whereas 
Caprion uses a Lorentzian fit. It was noted earlier that the peaks in S(q, uj) are asymmetric 
34j . This makes a fit dependent on details of the adopted fitting procedure. 



B. Structure factor at high (/-values 

The dynamic structure factor in the Brillouin regime provides only rather limited infor- 
mation on the structure of the vibrational modes, even if one includes the experimentally 
hardly measurable transverse branch. Additional information can be gained if one includes 
q- values beyond the first diffraction peak where S(q, uj) can be measured by inelastic neutron 
scattering. As seen from Fig. 0] the dynamic structure factor exhibits even for this simple 
material a rich structure. For the sake of clarity we plotted Sft((?, / q 2 Z {uj) . Taking con- 
stant uj cuts of S{q,uj) one gets distinct patterns for the average of the modes with that 
frequency. 

In Fig. we show four such constant omega scans: for uj = o>bp, for the top of the 
apparent acoustic phonon in Fig. EJ for uj near the Einstein average and for localized vibra- 
tions. For comparison the static structure factor, S stat {q), is also shown. Comparing the low 
frequency and high frequency vibrations one observes a clear shift of phase with respect to 
the static structure factor. The low frequency vibrations are in phase but the localized ones 
are in anti-phase to S sta _ t {q). In simple densely packed materials, e.g. metals, the nearest 
neighbour coupling dominates. Localized vibrations are vibrations of neighbouring atoms 
against each other and this is seen as phase shift against S' stat (g). Low frequency vibrations 
on the other hand must not involve significant changes of the nearest neighbour distance 
and therefore show no phase shift. The intermediate frequency modes {uj =3.8 and 8) show 
a mixed behaviour. The lower of the two frequencies corresponds to the high frequency 
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FIG. 4: Dynamic structure factor Sft(q, uj)/q 2 Z{uj). For oj = the static structure factor was 
taken. 

edge of the transverse mode. In the constant-omega scan one sees a narrow line at low q 
from the longitudinal phonon, a flat maximum around ofp and a weak modulation, clearly 
distinct from S' stat (g) at larger q. For the other scan with oj =8, which is above the apparent 
maximal transverse phonon frequencies, one sees again the narrow peak of the longitudinal 
phonon, now at a higher q- value, and additionally a broader peak marking the descent of the 
pseudo-dispersion for q — > qFP. The oscillations for larger q are out of phase with S sta t{q) 
and damped compared to the curve for the localized vibrations. Clearly the character of 
the modes with oj « c^bp is distinct from the one of the zone boundary modes. Using our 
previous results, we take a closer look at the modes responsible for the BP. 

In Ref. (34 1 we deconvoluted the exact low frequency eigenvectors of the present system 
into extended sound wave like modes and local modes with a weak bilinear interaction. 
This was done by rotating the basis of modes in a narrow energy range. This way we 
extracted optimally localized vibrational modes and automatically gained the proper number 
of extended phonons. The total number of modes is preserved and the number of modes in 
a given frequency interval is only marginally altered. The resulting extended and localized 
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FIG. 5: Constant frequency scans of the dynamic structure factor. Solid line: static structure 
factor, a: dashed line u = 0.6 f» lubp, long dashed line to = 3.8, dotted line uj = 8, dash dotted line 
w = 12; b: uj = 0.6 w wbp, dashed line averaged as in a, dotted line: extended mode, dash dotted 
line: local mode. 

modes, are no longer eigenvectors of the dynamic matrix but have a weak bilinear interaction. 
This is the picture used in the soft potential model. In Fig. EJd we compare S sta x{q) with 
constant-u; scans at ~ c^bp for the the true harmonic eigenmodes and the deconvoluted 
extended and local parts. It is not surprising that the curves of the deconvoluted extended 
modes follows 5' stat (g) even more closely than the one of the eigenmodes. This property 
of sound was was suggested by Buchenau et al. |3] as a distinguishing feature between 
them and local vibrations. This difference is clearly not very pronounced in the present 
case. For the deconvoluted local low frequency modes the first peak is certainly reduced but 
still quite pronounced. For larger q the scans of the extended and local modes are hardly 
distinguishable. This somewhat surprising result is a consequence of the strongly collective 
nature of the local vibrations, which form the cores of the quasi- localized vibrations seen 



in Excess- These modes are collective motions of chains of 10 and more atoms 
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Such modes are typical fore densely packed metallic glasses. 



IV. RELAXATIONS 

Besides vibration one observes in a glass aperiodic jumps of groups of atoms, local relax- 
ations. Such jumps give e. g. rise to the t eleg raph noise in point contacts and are thought 
to be the elementary process in diffusion 50]. In our previous study of relaxations in the 
soft sphere glass we found them to be closely correlated to the quasi- localized vibrations. 
They are again collective motions of chains of atoms (44, 49]. At low temperatures any 
given atom moves in such a jump only a fraction of the nearest neighbour distance. One can 
observe these collective jumps in a molecular dynamics simulation by monitoring the total 
displacement of all atoms as function of time 



AR(t)= Y,(mt)-mo)f (9) 

y n 

where R n (t) is the position vector of particle, n, at time, t, and R n (0) is the one at the 
starting or reference configuration. AR(t) oscillates due to the vibrations and changes due 
to relaxations, i.e. due to the transitions from one local energy minimum to another. At low 
temperatures T pa 0.05T g , the maximal distance an individual atom travels in these jumps 
is only 0.3a, about a quarter of Rnn- Increasing the temperature by a factor of four the 
jumps seen at the lower temperature can no longer be resolved and new jumps over larger 
distances are observed. 

All observed relaxations were collective jumps localized to 10 and more atoms forming 
twisted chain-like structures with some side branching 44]. This collectivity is again reflected 
in the relaxational structure factor, defined in analogy with the vibrational (Eq. Ej) one 



S'rel(Q) - 



^QAR"exp(?QR r ' 



(10) 



where AR" is the jump vector of atom n and (. . .) denotes averaging over angle and re- 
laxations. The relaxational structure factor is similar to the corresponding curve for the 
low frequency vibrations. It follows its static counterpart closely, independent of the jump 
lengths. The smoothing out of the higher peaks is more pronounced for the longer jumps 
occurring predominantly at higher temperatures. Comparing Figs. and |B1 one sees the 
close correlation between quasi-localized vibrations and atomic jumps expressed which we 



established earlier by correlating vibrational eigenvectors an jump vectors 
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FIG. 6: Structure factor of local relaxations at temperatures up to 0.15T g . Full line: static 
structure factor, dashed line: average over relaxations with AR < 5a, dash dotted line: average 
over relaxations with AR < 10a. 

V. CONCLUSIONS 

The vibrational spectrum of the l/r 6 -soft-sphere glass shows a pronounced boson peak. 
The boson peak frequency is very low at about 5% of the maximal frequency. The shape 
of the boson peak is in accordance with the one obtained from a model of interacting soft 
oscillators and sound waves. Calculating pseudo dispersion curves one finds a different 
behaviour for longitudinal and transverse phonons. Whereas the longitudinal ones are well 
developed in the whole of the first pseudo Brillouin zone the transverse ones become over- 
damped near the boson peak frequency. Perhaps not surprisingly, but often forgotten, this 
shows that depending on the underlying atomic structure the two branches of sound waves 
vary in their behaviour as function of q. In the studied model glass, the boson peak vibrations 
have predominantly transverse character. 

From the structure factor at higher q we find that the boson peak vibrations are distinct 
from the zone boundary modes. They closely follow the structure of the static structure 
factor. Splitting the exact eigenmodes into extended phonons and localized modes one finds 
that both components again follow the large q part of the static structure factor. The larger 
extent of the phonons is seen as a larger amplitude at the first peak. This semblance of 
the structure factors can be understood from the collectivity of the localized vibrations, 
the cores of the quasi-localized vibration. As we found earlier these modes are collective 
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vibrations of chains of ten and more atoms. These chains move in such way that the nearest 
neighbour distances do not change markedly. 

The structure factors of local relaxations again resembles the static one and is closely 
related to the one of the boson peak vibrations. This is seen as a consequence of of the 
similar structure of the low frequency localized vibrations which produce the boson peak 
in inelastic scattering and the collective jump processes which are the elementary step in 
diffusion and relaxation in metallic glasses. 
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